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Abstract. We consider the phase diagram of twisted mass Wilson fermions of two- 
flavor QCD in the parameter space of the quark mass, the isospin chemical potential, 
the twist angle and the lattice spacing. This work extends earlier studies in the 
continuum and those at zero chemical potential. We evaluate the phase diagram as 
well as the spectrum of the (pseudo-)Goldstone bosons using the chiral Lagrangian 
for twisted mass Wilson fermions at non-zero isospin chemical potential. The phases 
are obtained from a mean field analysis. At zero twist angle we find that already an 
infinitesimal isospin chemical potential destroys the Aoki phase. The reason is that 
in this phase we have massless Goldstone bosons with a non-zero isospin charge. At 
finite twist angle only two different phases are present, one phase which is continuously 
connected to the Bose condensed phase at non-zero chemical potential and another 
phase which is continuously connected to the normal phase. For either zero or maximal 
twist the phase diagram is more complicated as the saddle point equations allow for 
more solutions. 


1. Introduction 

In the past decade lattice QCD simulations with twisted mass Wilson fermions have 
attracted a great deal of attention [T]-[^. The advancement of algorithms as well as the 
increasing power of computers have allowed for simulations in the deep chiral regime. 
Simulations with ordinary Wilson fermions in this regime can be severely hindered by 
the so-called exceptional conhgurations. These are conhgurations for which a Dirac 
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eigenvalue is extremely close to minus the quark mass so that the Dirac operator cannot 
be inverted. With a twisted mass on the contrary the determinant is regulated by 


the mass, see for example (1.3) below. The twisted mass is a mass which comes with a 

for excellent reviews of the twisted 




75 in Dirac space and a ts in flavor space (see 
mass formulation). In the twisted mass basis the fermionic part of the Lagrangian reads 

.p =+ zu'y5T3)'ip = + m + zu'y5T3)'ip, (1.1) 


where Dw is the Wilson Dirac operator, m the quark mass ( i.e D^{m) = Dw + the 
massive Wilson Dirac operator), and u the twisted quark mass. The quark fields are 
denoted by x and y. One can perform the following change of variables x = 
and X = where one introduces u = arctan(M/m) and M. = \/m^ + u^. This 

transformation defines a new basis, the physical basis, where the Dirac operator assumes 
the usual form 


4” = x(^w + -^)x. (1.2) 

However as indicated by the prime, the Wilson term is now rotated in the Dirac operator. 
Since the Wilson term vanishes in the naive continuum limit, twisted mass QCD should 
be equivalent to ordinary QCD. The proper mapping between twisted mass QCD and 
ordinary QCD was analyzed in . The addition of the twisted mass renders the fermion 
determinant positive definite. 


det{Dy/{m) + iuj^ts) = detdet{Dy/{m) — iu'^3) 

= det(Dw(m) + zttys) det(75Dw(m)75 - lu'y^) 

= det(Dw('m)D^(m) + M^), (1.3) 

where we used the 75 -Hermiticity of Dw = The other crucial property of 

twisted mass Wilson fermions is the automatic 0{a) improvement of both the actions as 
well as of the operators which is achieved at maximal twist without significant additional 
computational cost compared to standard Wilson fermions [^ . It is noteworthy to point 
out that the case of intermediate twist is not of pure academic interest. In [9,37,38 


intermediate twist angles were used to achieve 0{a) improvement in the simulations 
dedicated to the generation of ensembles employed in renormalization projects. The 
simulations with very small PCAC fermion masses were numerically unstable and had 
very large autocorrelation times. This could be cured by averaging the results for 
opposite PCAC fermion masses over the twist angle, which also achieved an 0{a) 
improvement. This method was proposed in [^. 

In order to make full use of the results of lattice simulations with Wilson fermions it 
is essential to understand the non-trivial phase diagram which is caused by the specific 
choice of discretization. The phase diagram of Wilson fermions with a non-zero twisted 
mass has been studied in 10 -26 . In the present work we extend these studies to include 


the isospin chemical potential. Thus we generalize finite density studies in continuum 
QCD and QCD like theories 27 35 . QCD at finite isospin chemical potential is of 


interest in different physical contexts. Two prominent examples are neutron stars and 
heavy ion collisions. Moreover, studying QCD at finite isospin chemical potential can 
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give useful insights into the notorious sign problem in QCD at non-zero baryon chemical 
potential (^ . 

In section we recall some properties of the low energy effective theory, especially 
the p-regime, for twisted mass Wilson fermions. Thereby we present the three leading 
order terms of the chiral Lagrangian and introduce the order parameters. The 
thermodynamic limit is performed in section We start the discussion of the phase 
diagram with an analysis of the case with vanishing twist {u = 0) and determine the 
different phases that one encounters in the cases of zero, real and imaginary isospin 
chemical potential. Then the phase diagram can be studied in the plane of the quark 
mass and the lattice spacing. At zero twist we find that the Aoki phase only occurs 
at zero isospin chemical potential. Thereafter we also discuss the structure of the 
phase diagram at finite twist. Finally, we consider the phase diagram of the physically 
important case of maximum twist. The detailed computations to this analysis are carried 
out in [Appendix A 


In section]^ we discuss the dependence of the pion masses on the chemical potential 
and lattice spacing. The detailed calculation to this part are given in Appendix B In 
section |5] we summarize our results. 

The present work was partly published in the proceedings [^. Moreover we want 
to point out that we work with natural units (c = h = /cb = !)• 


2. Chiral Lagrangian of Twisted Wilson Fermions 

We consider the partition function of two dynamical Wilson fermions in the fundamental 
representation of the gauge group SU(3) with degenerate quark masses m. 


/ D|C/]exp|-Sg,„j,(C/)] det (£)„( 7 n) + y7oT3^ , 


( 2 . 1 ) 


where 5'gauge(f^) is a discretized version of the Yang-Mills action, e.g. the Wilson gauge 


action. Here we follow the notation in 29,39 for the isospin chemical potential which 


enters with a factor 1/2. We consider both real and purely imaginary isospin chemical 
potential /ij. 

We emphasize that the ensuing discussion also applies to SU(Ai/ > 3) with the 
fermions in the fundamental representation. Those theories have the same pattern 
of spontaneous symmetry breaking as the one of QCD, and hence, share the same 
chiral Lagrangian for the Goldstone bosons, though the low energy constants might be 
different. 

The operator is the Wilson Dirac operator with a twisted mass, defined in 


(2.4) below. Its properties are briefly recalled in subsection 2.1 In subsection 2.2 


we introduce the chiral Lagrangian in the p-regime as well as its mean field limit 
which coincides with the leading order Lagrangian in the e-regime also known as the 

The order parameters which become important later on are 
They are given in terms of the parametrization of the 


microscopic limit 40 


presented in subsection 2.3 


Goldstone manifold SUl(2) x SUr(2)/SUv(2) = SU(2) in subsection 2.4 
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2.1. Properties of the Euclidean Wilson Dirac operator with twisted mass and isospin 
chemical potential 


Let 7 q, be the Euclidean Dirac matrices. Then the Euclidean Wilson-Dirac operator , 

DY/{m) = 7'^Dk — aV^ + m, (2.2) 

is 75 -Hermitian, = 75 Dw 75 - The Dirac operator with isospin chemical potential 
is 75 = diag (1,1, —1, —1) Hermitian if /xi is imaginary otherwise it becomes 75 ri /2 
Hermitian, i.e. 

75 D (^Dw{m) + y7oD) 75D = 75 D 75D, 

/ T-. / \ hi \ ^ hi ^ 

(^Dw(m) + Y7oDj = 


75 75, 


Pi G 


( 2 . 


From both Hermiticities it follows that the spectrum of D\^{m) + (pi/ 2 ) 7 or 3 consists of 
complex conjugate pairs or real eigenvalues. However the product of the real eigenvalues 
is positive only if pi is real since {'^ 5 D^^{rn) + (^ 1 / 2 ) 7570 )'*' = 'y 5 Dy/{m) — (^ 1 / 2)7570 so 
that A*(pi) = A(—Pi). This guarantees the positive semi-dehniteness of the statistical 
weight in Eq. (2.1) if pi G M while for imaginary pi only the reality of the statistical 


weight is ensured. 


The twisted mass QCD partition function (2.1) with the Dirac operator is given by 
D^{m) = — aw + mcosu + imsmu'-y^r^. (2.4) 


Any kind of Hermiticity of D^ is lost for an arbitrary twist oj if the isospin chemical 
potential is imaginary. Thus the positive dehniteness of the statistical weight will never 
be achieved in this case. However, the Dirac operator is still 75 ri /2 Hermitian if pi is 
real such that the statistical weight is still positive semi-dehnite. 

In the ensuing discussions we restrict the twist to a; G [0,7r/2]. Indeed the twists in 
the other three intervals [—tt, — 7r/2], [—7r/2,0], and [7r/2,7r] can be traced back to the 
case u G [0,7r/2] by relabelling the up and down quarks u ^ d and/or reflecting the 
quark mass m —)■ —m. Therefore this assumption does not restrict generality but it has 
the advantage that both cos a; and sin a; are non-negative. ‘ 


2.2. Chiral perturbation theory of Twisted Mass Fermions 

The chiral Lagrangian of Wilson fermions with a twisted mass and a hnite chemical 
potential follows from the symmetries of the fermionic action. Classically it is 
invariant under the group SUL(A^f) x SUr(W) x Uv(1) x Ua( 1)- Due to the chiral 
anomaly the group Ua( 1) is broken explicitly. The chiral symmetry is spontaneously 
broken by a non-zero chiral condensate S = |(V’' 0 )l with symmetry breaking pattern 
SUL(fVf) X SUR(iVf) —t SUv(Af). Therefore the corresponding Goldstone bosons U are 
elements in f/ = f/oexp[zn/Yr] G SU(iVf) with H the pion helds. 






Phase Diagram of Twisted Mass Wilson Fermions 


5 


To construct an effective theory one needs two ingredients. The hrst ingredient 
is the thermodynamic and low-mass, low-momentum, low-chemical potential and low- 
temperature limit so that the QCD partition function is dominated by the pseudo- 
Goldstone bosons. Second, we need a counting scheme to order the terms in the chiral 
Lagrangian allowed by symmetry. We choose the p-counting scheme, 

~ 1/1/, /if ~ 1/G, ~ 1/1/, pf ~ 1/1/, ~ 1/1/ hxed, (2.5) 


where a is the lattice spacing and p^ is the four momentum of the pion helds 11. In this 
counting scheme the kinetic modes (p 7 ^ 0 ) and the modes with zero momentum (p = 0 ) 


are still coupled. To lowest order the action is given by 11-15 


5p = / d^x Cp{U{x)) 


( 2 , 6 ) 


with the chiral Lagrangian 

CpiU) = ^tvD^UD^U^- 


mT, cos a; 


tT{U + U^) 


imT sin oj 


trr3(17-t/t) 




The covariant derivatives in flavor space are given by 27,28 


D,U = d,u - ^[f/,r3]_<5,0, 


hi. 


(2.7) 


( 2 . 8 ) 


with [.,.]- the commutator. The symmetries also allow a term linear in the lattice 
spacing a. However, this term is proportional to tr(f7-|- U^), such that one can eliminate 
it by renormalizing the quark mass m and the twist angle 00 . Therefore we can omit 
this term without loss of generality. 

The low energy constants are the chiral condensate S, the pion decay constant 
Fjr and the two low energy constants hhe/s whose sign convention is chosen as in 
Refs. 42-45 . The third low energy constant associated with the of terms, usually 


denoted by — f/f) vanishes for the two flavor SU(2) partition function, see 

Ref. 20 . However this does not mean that this low energy constant does not play 


any role in two-flavor QCD. It still affects the eigenvalue spectrum of the Wilson Dirac 

for analytical discussions of this operator at vanishing isospin 


43 45 


operator, see 
chemical potential. 

The two low energy constants corresponding to the hnite lattice spacing reduce 
to one constant C 2 = Wq -|- H 8/2 because the trace and the squared trace of a matrix 
U G SU(2) are related to each other, tr^ U = tr f/^ -|-2. While W% is positive dehnite, see 
the discussions in 26,42,46 , the sign of Wq is negative so that C 2 can have either 
sign. 

We will derive the phase diagram using this effective theory. In order to do so it is 
useful to introduce four additional source terms which probe the vacuum structure 
3 

^ ^ jk^lsTk'tp + 

k=l 


(2.9) 
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which in the chiral Lagrangian lead to the terms 


C,onrce{U) = _ f/t) _ tr(17 + U^). 


( 2 . 10 ) 


Note that the jo"term and the mv-terms are not really new source terms since they 
enter like the twisted mass terms mcoso;! and ^msincuysTs. Thus they can be 
certainly absorbed by the quark mass m and the twist angle u. Nonetheless we still 
include these two terms due to two reasons. First, we want to underline the SO(4) 
transformation property (the adjoint action of the SUl(2) x SUr( 2) flavor group with 
tjj —)■ exp[?(Q;^rfc+ /3^75rfc)]'^ and tjj —)■ exp[?(—+ /3^75rfc)] where and (3^ are real 

angles) of the four dimensional Euclidean real vector j = {ji, j 2 , Second, the 

sources jo and can be chosen independently of the quark mass m and the twist angle 
u. This becomes particularly important when either ca = 0,7r/2 or m = 0 and the saddle 
point manifold is more than only a single element of SU(2) in the thermodynamical limit. 

The phase diagram will be determined from the saddle point where the effective 
action Sp takes its minimum. Let Uq be this saddle point. Then the held Uq has to be 
constant over the whole four-dimensional box V since each kinematic part will increase 
the action. Therefore we expand the held U as 


U = Uoexp A^n(p,E) 


exp {i{pkX^ — Ef)) 


( 2 . 11 ) 


= Uo ll 2 + ^ 5 ^n(p,E) 


exp [i{pkX^ — Et)) 


W n(p. EW, ~ + o(r-^/^) 

f ^ 1 


V.E.'p'.E' 




-%[pkX^ 

dTT^ 


( 2 . 12 ) 


5: n.(p. t,. + o(r-7. 


^ p,E,p',E' J 

where we employ the Fourier transform of the pion helds. Note that the unitarity of U 
enforces the relation 

U{p,E) = Il\-p,-E). (2.13) 

The factors of l/27r are the normalization factors of the plane waves such that we have 


VS(E-E')llS(pt-A)- 


(2.14) 


To improve the notation we dehne the dimensionless quantities 


a2 = aVc2, f = p^VEl E^ = E^VE^ 


in = mPV, pI = p^VE^, 


( 2 . 15 ) 
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Then the expansion (2.11) is plugged into the chiral Lagrangian (2.7). Expanding up 


to second order in 11 the action reads 

j4 


= vCoiUo) + j d^x £i(t/o,n) + y d^x £2(t/o,n) 

+ j d^X £source(t/o) + (2.16) 

where £ 0 (^ 0 ) is the mean field chiral Lagrangian evaluated at the saddle point Uq 


VCoiUo) = ^ tr[f/o, TsUUl ^ 3 ]- - Y tr(e*""3^o + tUiUo + U^o)- (2-17) 


To the same order we have to evaluate the source terms C 
terms CoiUo 
by 


fUo)- The zeroth order 


£source(7^o) are of order 0(1/ while the first order terms in 11 given 


d^a:£i(f/o,n) = 


^ fn( 0 ) f (e^-’C/o - C/Je-””) 


+ 23 "(f/„-f/„')tr(f/„ + ( 7 |J 
and the second order terms in 11 given by 
2 


( 2 . 18 ) 


j A A(c/„,n) = /j ^ 


p,E >0 L 




£2 ^ ^ TlfE ^ -^4. 

— trn(p,E)nHp,E) - ^triUon(p,E),nHp,E)u^].T3 


+ ^Re tr[n(p, E), r3]_[nt(p, E), U^TsUo]. + ^ tr(U^^^Uo + Uy^^^=^)n(p, E)nt(p, E) 

lo 4 


-\-(l 


tr(f/„ - r/„')n(p, E) - S'" tr(f/„ + c/t[)n(p, B)nt(p, e) ti{u„ + u^) 


( 2 . 19 ) 


are of order 0(1/“^/^) and 0(1/“^), respectively. 

The saddle point Oq is calculated in [Appendix A| and will be discussed together with 
the corresponding phases in section At the saddle point Uo the first order Lagrangian 
Cl vanishes. The Lagrangian C 2 is the lowest non-vanishing order for the non-zero 
momentum modes of the pseudo-scalar pions. In section we will derive the dispersion 
relations using £ 2 ( 00 , If) and subsequently read off the masses. This is also the reason 
why we Wick rotated from Euclidean to Minkowski space, i.e. do —)• —idt. 

Before we proceed with this discussion, let us note that the shifted Lagrangian 


VCoiUo) — 8W, see eq. (2.17), is invariant under the following transformation 


TT 

oj —y — — CJ, 

2 

Uq TiUoT2^ 
ct —y — ct , 

/Xj — y /Xj -|- 32(X . 

To see this one has to use that any Uo € SU(2) can be parametrized by 
Uo = al 2 + tfdkTk with + (dkPk = 1 and a, ftk G M. 


( 2 . 20 ) 


( 2 . 21 ) 
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Since this symmetry of £o is valid for all Uq G SU(2) it can be extended to £i with the 
additional change n( 0 ) — r 2 n( 0 )r 2 because Ci is a derivative of Cq. The symmetry is 
violated by the kinetic term (in particular the term proportional to fliE) and does not 
hold for £ 2 - Therefore the pion masses do not reflect this symmetry. As we will see 
below, this symmetry relates the phase diagram for maximal twist to the phase diagram 
for zero twist. 


2.3. Order Parameters 


The phases of the full partition function are determined by 

Zo = [ dniUo) exp[-V{Co{Uo) + Csonrce{Uo))] 


( 2 . 22 ) 


and are characterized by order parameters. We consider the following ones: 

i) The chiral condensate can be introduced by introducing the auxiliary variable m^, 

^ 11 /'- 
S(m) = — - lim — / dl^x{'il){x)il){x)) 

2 ,v^ooV J 

= lim lim — 77 ^—InZo 
mv-5>0 V-s-oo 2V OUlv 


= 7Jim {tT{Uo + Ul)). 
4 


(2.23) 


Note that only is a renormalization group invariant, and the actual order 

parameter is given by this combination. In the case of the Sharpe-Singleton scenario, 
the quark mass has to be replaced by the distance between the quark mass and the Dirac 
spectrum. Since it remains finite for m —)■ 0, we obtain a first order phase transition. 
Strictly speaking it is not a phase transition because the two phases have the same 
physical properties. In particular the system can lie in two different but physically 
equivalent ground states at vanishing quark mass, 

ii) The 7 r° condensate can be generated by the source jo = jo, 

i If 

= - lim — / d‘^x{%lj{x)'-^oTo'ilj{x)) 

Z V^oo y I 


= lim lim —— t^IuZq 
jo-s>OV^>oo 2V ojo 

= l^{tlTo{Uo - ul)). 


(2.24) 


iii) For the charged pion condensates we employ the standard notation that oc du and 
n~ oc ud such that tf = {u, d). Therefore when dehning = Ti ± 1 x 2 and j± = ji ± ij 2 
these condensates read 

= lim lim —— 7 -—IuZq 

>-oo 2V oj± 


(2.25) 
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iv) The isospin charge density is given by 


ni 


V-5>oo ZV I 


InZn 


lim 

V^OO 2V Ofil 

-^(trlC/o.Tsl-lC'o'.Ta]-). 


(2.26) 


Since the phases are already fully characterized by the hrst three condensates we do not 
need the charge density as an additional order parameter. 

Nonetheless other interesting observables such as the pion covariances and the chiral 
variance always play an important role when the angles of the Goldstone manifold SU(2) 
do not completely freeze out but the ground states are only hxed via the source terms. 
These three covariances are given by 


AG^o = lim 


1 




= — lim lim 


jo-s>o y-5>oo ( 9 j| 


A ~ 

X1p{x)'^^T‘j,f){x] 

rlnZn 


d'^X'lp{x)-f5T3'lp{x] 


16 L 


- U^o)) - (trr3(17o - Gj)) 


t ^\2 


(2.27) 


for the variance of the 7r° condensate, 


AG.,r± = lim 


y^oo 41/2 


= lim lim 






dh^xil){x)'^3—il){x 




i±-5>o y^>oo 41/2 

s2 r /i , i2 


InZn 


64 


tiT+iUo-Ul) 


trr+(f/o - Ul) 


for the covariance between the two charged pion condensates and 

f - ^\ / f - ■' ^ 

dfxil){x)il){x) ) — { dfxil){x)il){x\ 


AS(m) = lim 


y^oo 41/2 
1 


= lim 




y^>oo 4 G 2 Qrnl 


16 L 


InZo 

lmv=mcoscj 

rt\\ _ _L rrt'i \2 


(tr2([/o + GJ)) - (tr(f/o + Gj))^ 


(2.28) 


(2.29) 


for the variance of the chiral condensate. We want to emphasize that these covariances 
would be proportional to the susceptibilities at hnite volume V, i.e. our dehnition 
incorporates an additional factor of 1/1/. However in the case that the average 
of a condensate vanishes the corresponding (co-)variance remains hnite while its 
susceptibility is of the order OiV). 



Phase Diagram of Twisted Mass Wilson Fermions 


10 


2.4- Parameterization of Uq 


The mean field phase diagram is determined by the zero momentum part of the chiral 
Lagrangian (2.17) with parameters occurring in the combination 

(2.30) 

The terms of the zero momentum chiral Lagrangian coincide with the leading order 
chiral Lagrangian of the e-counting scheme as well as with the chiral Lagrangian that 
is obtained from chiral random matrix theory in the e-domain of QCD. Therefore the 
mean field phase diagram can also be obtained from the Wilson random matrix partition 
function 42 . However, to determine the pion masses we need the kinetic term and have 
to use the full lowest order chiral Lagrangian in the p-counting scheme. 

To perform a saddle point analysis we need a parameterization of the group SU(2). 
Our choice is 


Un = V<!>V 


-1 


(2.31) 


r 

cosy 

sin 

^ ■ 

2 



0 


e 2 sm — 

2 

cos 2 



0 






cos ^ 

-e ^sm — 
2 


sin 




cos 2 


= COS</9]l2 + ^sin<y9 


■ g^l?2/2 

0 


cos'll sin'di 


0 

£-*>52/2 


sin 'di — cos i9i 



0 


0 


with the angles (p G [0,7r], 791 G [0,7r[, and '^2 G [0,27r[. The Haar measure dp{U) in 
terms of these coordinates follows from the invariant length element, 

tidUdU-^= iid^d^-^+ ii[^,V-^dV]_[^-\V-^dV]_ (2.32) 

= 2dip^ -I- 8 sin^ ipd'd\ -|- 8 sin^ (p sin^ -did-dl, 
such that we have the invariant measure 

dfi{U) = 16sin^ 93sin'did'di(i'd2dv3. (2.33) 

Note that sinip as well as sin'll are positive on their support. 

In the coordinates (2.31) the chiral Lagrangian (2.17) is given by 

VCq{(p, "di) = —2m cos a; cos </9-|-2m sin a; sin ip cos'di-|-16a^ cos^ ^~~2~ ^ 

and 


Source( t^, "di) = S(2josin(pcos'll -f j+sinipsin'dic*’^^ -|-j_ sinipsin'dic *’^^) 

— 2EmyCos(p. (2.35) 

Note that the dependence on the angle '^2 completely drops out of the microscopic 
Lagrangian £o which determines the phase diagram. The vr^ condensates have a phase 
associated with an exact Goldstone mode corresponding to the angle i92- The exact 
nature of this mode can be fixed by the introduction of an appropriate source term 
which are the two terms proportional to j±. 
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In the thermodynamic limit the angles (p and "di are hxed by the saddle point 
equations. To measure these angles and, thus, in which phase the system is, we use the 
order parameters discussed in section 2.3 In terms of the parameterization (2.31) they 
are given by 

S(rh) = S lim lim (cos<p), 

m-v—>0 V^oo 

= S lim lim (sin(pcos'd!), 

jo—^■O V—^■oo 


= -lim lim (sincpsin'dic^*’^^). 


(2.36) 


Note that the order of the derivatives and the thermodynamic limit is crucial since they 
do not commute. The isospin charge density in terms of these variables is given by 

WFI 


ni = :-^ lim (sin^ (psim'di). 

2 n^oo 

Also the variances take simple forms like 

lim ((cos^'di sin^ (p) — (cos'd! sin(p)^) 

y^oo ^ ' 

for the 7r° variance and 


(2.37) 


(2.38) 


AE = T? lim ((cos^(p) — (cos(p)^) 

y^-oo ^ ' 


(2.39) 


for the variance of the chiral condensate. The covariance between the two charged pion 
condensates can be traced back to the isospin density, 

mTf 


Aa± = 




la 


(2.40) 


and, hence, does not yield anything new. 

In the case that the angles <p, 'd! and 'd 2 freeze out at only one saddle point Uq due 
to the large volume limit 17 —)■ cxd (including source terms j ^ 0) all quantities reduce 
to two terms only, namely 


E(m) = Ecos<po and C^o = S sin (po cos'd!,o, (2-41) 

where the subscript 0 indicates the value at the saddle point. The other terms take the 
values 



AC^o = AC^i = AS = 0. 


(2.42) 


Note that the variances and covariances only vanish because the angles take certain 
values when switching off the source terms they can be non-zero. 
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To determine the phase diagram we have to minimize the effective potential given by 


VCq in Eq. (2.34). It is useful to rewrite the potential in the form 

2 


VCo = 


i6a= + 1 


2 m cos u 


cos (p 


2 


32a2 + 
2 m^ sin^ oj 2w? cos^ u 


^ 2 


sin ip cos'll + 


2 m sin oj 


(3.1) 


fp 32cP + pj ’ 

and we have to determine its minima with respect to the variables ip G [ 0 , vr] and 
■di G [0, 7 r[. The third angle '^2 parametrizing the group SU(2) can only be fixed by the 
source terms j±. 

From Eq. (3.1) one can easily read off two cases determined by the solutions of the 
saddle point equations for ■di, 

sin (p cos'd! = 


2 m sin u 
Ti 


and sin'di = 0 . 


(3.2) 


i) The hrst solution in Eq. (3.2) is assumed when the isospin chemical potential is real 
{jl\ > 0 ) and the modulus of the mass is bounded from above as 0 < 2 |m|sinC(; < 
pf sinip. The Lagrangian simplihes to 


VC.= 




cos ip — 


2m cos OJ 
32a? + 


2 


T 2rh? sin^ oj 2rh? cos^ u 


Ft 


32a? + 


Minimizing with respect to (p we hnd two different solutions 

2 m cos u 

cosip = ^ or sinip = 0 . 


(3.3) 


(3.4) 


3222 + jq 

The phase with 0 < cosip < 1 and 0 < cos'd! < 1 will be denoted by I. The 
solution sin(y 9 = 0 can only exist for msincj = 0 (denoted by 1//^'^^°^). At zero twist we 
also have the solution cos'd! = 0 for fi\ > 0. This gives two possible phases, one with 
I cos</ 9 | = 1 and one with | cos^?! < 1. Otherwise the hrst solution in Eq. (3.2) cannot 
be satished. For oj = 'n 12 we hnd cosip = 0 and cos'd! = —2m//if which will be referred 
to as phase / 0 =’r/ 2 )_ 

Further analysis of this case including the computation of the parameter domain 
of the various phases is carried out in [Appendix A.l 

ii) The second solution in Eq. (3.2) minimizes the ehective potential when 
2|m|sina; > pifsimp which includes the case pi < 0. Then cos'd! has to take one 
of its extremal values ±1. Comparing both solutions we hnd that cos'd! = —signm is 
always the minimum in these cases. Then the Lagrangian reads 

VCo = 16a^ cos^ ip — 2|m| sincasin^? — 2 mcosa;cos 93 . (3.5) 


This Lagrangian satishes up to a constant shift the following symmetry 

u ^ — u, if ^ —sign(m) — cp, a ^ —a , (3.6) 

which is reminiscent of the symmetry ( |2.20[ ). Note that no chemical potential is involved 
in this symmetry since the free energy is independent of it in this part of the phase 
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Figure 1. Solution of the transcendental equation (3.7) for ip S [0,7r/2], see Eq. (A.11) 
as an explicit function of |m|/16a^. Note that the solution is not differentiable for 
uj = 0,7r/2 which is the reason for the splitting of the phase II into two phases at zero 
and maximal twist. 


diagram and in this order of the free energy. Therefore the system enters the “Silver- 
blaze-property” . 

The saddle point eqnation can be rewritten as 

8 a^ sin 2(p = —\m\ sin(a; — sign {m)ip) (3.7) 

and can be solved nnmerically or analytically as we do in [Appendix A.2 In Fig. [T]we 
show the solntion for a range of twist angles. For 0 <ti;<7r/2 we always have that 
I cos</9| < 1 for I cos'll = 1. This phase will be denoted by II. 


For a; = 0 the solution of the transcendental equation (3.7) splits into two branches 


cos(p = 1 (denoted by phase and cos(p = m/16a^ (denoted by II^‘^^^'>). We 

underline that the angle i?! is not determined in the hrst case because the saddle point 
is Uo = 12- Also for u = 7 t/2 the solution branches into cosi^i = 1 with coscf = 0 
(denoted by phase IIlP='^D)'^ and cos'll = 1 with sin^ = m/16a^ (denoted by phase 

J-J-(ui=TT/2)y 

We further analyze case ii) 


m 


Appendix A.2 


We hnd two phases at hnite twist oj ^ ^■,t^I2. and three at each extremal twist 
uj = 0,7r/2. Additionally there is an Aoki phase at a; = 0 which has a counterpart 
at a; = 7r/2 due to the symmetry (2.20). The parameter domain of the phases is 


summarized in Table The corresponding order parameters are collectively presented 
in Table Thereby we also show the covariances which are non-trivial for those 
condensates which are only aligned by the source terms j and m^. We underline that at 
j, ruv = 0 the condensates vanish while the covariances do not. However the condensates 
do not vanish if the thermodynamical limit is taken before the limit j, —>■ 0. This is 
a fundamental characteristic of spontaneous breaking of symmetries. 
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phase 

region 

I, u e [0, tt/2] 

Rj > 2 m sinca > 0 (if a; = 0 then yUj > 0) 

1 --2 2uT m costa ^ ^ 

nnrl -1- ^ ill ^ 9 m, ms//) t> 0 

\ \ 1 1 ti \ '.'V KJ 

y/fif - Am? sin^ ca 

II±, oj e]0,7r/2[ 

^0 ^ , 2 m uiP costa 

//? ^ 9 7D 5^9/7 -1- //, • 

IJjJ \ t ! L \ Olll Uv \J CtXll-l , - Kj^dKAj 1 /R'T 1 

“ v/R^-4m2sin"ca 

or 2 m sinta > Rj 

II 

yuf < 0 and 16a^ > m > 0 

IIlP=o} 

Rj > 0 and 2 m > 32a^ + flf; 
or 0 > Rj and m > 16a^ 

Aoki(‘^=°) 

/if = 0 and 16a^ > m > 0 


/if < 0 and — 16a^ > m > 0 
or —16a^ > m > 0 and — 32a^ > /if > 0 

j j 

/if < 0 and m > —16a^; 
or 2 m > yuf > 0 and m > —16a^ 

Aoki("^=^/2) 

/if = —32a^ and — 16a^ > m > 0 


Table 1. Regions of the phases. Note that the phase III± only occurs for w 7 ^ 0, 7 r/ 2 . 
The Aoki phase only exists at /ii = a; = 0 which connects the phases I to along 

the chemical potential axis via a hrst order phase transition. It has an analogue at 
imaginary effective lattice spacing (fp < 0) and maximal twist (w = 7 r/ 2 ) with real 
isospin chemical potential equal to hi = —32cP. This phase exhibits a similar behavior 
of the corresponding condensates. 


3.1. Zero Twist 

First we concentrate on the case of zero twist {u = 0). In this case the mean held 
Lagrangian is given by 

VCo\i^=o = —2mcosip + 16a^ cos^ V? — 4/i^ sin^ 99 sin^ di. (3.8) 

Before we discuss the general case of this kind of the Lagrangian we hrst discuss the 
cases of zero lattice spacing and zero chemical potential. 
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Table 2. 


List of order parameters in the corresponding phases. The function 


F(|TO|/16d^) used for the phases II± is defined in Eq. (A.IO). The subscript “±” 
refers to the sign of the quark mass rh. The charged pion condensates and the isospin 
charge density are C^± = -j±/{2\j±\)^W - T,^{fh) - C^o and m = 2giFlCl±/T?, 
respectively, cf. Eq. (2.42). They vanish when YP{m) + C^o = W. The quantities 
proportional to the sign of the source terms j and ruv would vanish if we would 
set j,mv = 0. Only due to the limit j, rUv —t 0, after the thermodynamical limit, 
those condensates acquire a non-vanishing value and reflect spontaneous breaking 
of the corresponding symmetries. This behavior is also reflected in non-vanishing 
vari ances. The one for the charged pion condensate is always = ni/{2giFf), see 

Eq. (2.40). The prefactor 1/2 in the variance of the phases results from 

averaging over the two minima while the factor 1/3 in the Aoki and Aoki-like phase 
results from averaging over the angle di (note the non-trivial measure (2.33) and the 
symmetry (2.20)). 
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3.1.1. Zero Twist and Zero Lattice Spacing. In this case we want to recall the phase 
diagram in the low-temperatnre limit as a fnnction of the isospin chemical potential and 
the qnark mass. The physics of this problem is well nnderstood, see 28,29 . The phase 
diagram which is shown in Fig. [^for of —)■ 0. At /ii = there is a phase transition to 
a Bose condensed phase of charged pions with one exactly massless Goldstone boson. 
For /ii < ttItt the masses of the psendo-Goldstone modes are given by 


m. 


= = 0 ) + qiai = 


/2Em 

~w 


+ QTh 


(3.9) 


where q is the isospin charge of the pions. 

In Enclidean space-time the interpretation of the chemical potential is that the 
retarded propagation is enhanced by exp{qfiiXo) for particles with positive charge q 
and snppressed by a factor exp(—g/iixo) for the anti-particles. Asymptotically, the 
propagator in the imaginary time direction (temperature) behaves as 

This behavior is determined by the pole mass which can also be obtained from 
the analytic continuation of the propagator to Minkowski space time. At non-zero 
temperature this results in a net non-zero particle density while at zero temperature, 
uninhibited propagation of the lightest pion only occurs for /ij > tUtt. 

At imaginary chemical potential in Euclidean space time, retarded propagation is 
modihed by a factor exp(ig/iiXo) and a factor exp(—zg/riXo) for advanced propagation. 
This corresponds to an imaginary pole mass. At imaginary chemical potential the 


partition function becomes Roberge-Weiss periodic 48 with period /xi = 27riT/3 where 
T is the temperature. In the zero temperature mean held limit (T = 0) we thus 
encounter this periodicity immediately. In the microscopic domain, though, we have 
/ii ~ l/^A^ which is well within the hrst Roberge-Weiss period |^. We will study the 
phase diagram in terms of the microscopic scaling variables, which have to be taken 
much larger than unity so that the saddle point approximation can be justihed. In this 
way we can consider an imaginary chemical potential while not running immediately in 
the Roberge-Weiss periodicity in the low temperature limit. 


3.1.2. Zero Twist and Zero Isospin Chemical Potential. At vanishing chemical 


potential the particular structure of the phase diagram was hrst predicted by Aoki 10 
and was studied in detail by Sharpe and Singleton 11 in the p-regime of chiral 
perturbation theory. 

For imaginary effective lattice spacing a the saddle point in p is at p = 0 or at 
(p = TT depending on the sign of the quark mass. The chiral condensate is independent 
of the mass but hips sign when the quark mass traverses zero. The two phases are 
physically identical but since the ehective potential is nonvanishing for m = 0, they 
are not continuously connected as a function of the quark mass. This is known as the 


Sharpe-Singleton scenario (ll 
subscript rehects the sign of t 


We denote these two phases by 
le quark mass (signm = ±1). 


where the 
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Figure 2. Phase diagram at vanishing twist (w = 0) in the m-o^-plane (Figs, (i), (ii) 
and (iii)) and in the /tj-m-plane (Figs, (iv) and (v)). The generic cases are: (i) zero 
isospin chemical potential (^j = 0), (ii) imaginary isospin chemical potential {gf < 0), 
(iii) real isospin chemical potential ()lj > 0), (iv) imaginary effective lattice spacing 
(a^ < 0), and (v) real effective lattice spacing (a^ > 0). The Aoki phase only exists 
at vanishing isospin chemical potential. Recall that tP is a combination of two low 
energy constants. Therefore we can have < 0 although the physical lattice spacing 
a is positive definite. The abbreviations S.S.^ and S.S.’^ refer to the original Sharpe- 
Singleton scenario where the chiral condensate jumps when the quark mass crosses 
the origin and a Sharpe-Singleton-like scenario where the condensate jumps. The 
subscripts ± corresponds to the sign of the quark mass. The tt^ condensate is only 
non-zero in the Aoki phase if a twisted mass (jo or a; 0) is introduced as a source 
term. This condensate spontaneously breaks parity and flavor symmetry. 


At real effective lattice spacing (a^ > 0) and vanishing chemical potential the chiral 
condensate does not jnmp anymore. When decreasing the qnark mass one enters the 
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Aoki phase through a second order phase transition. In this phase the saddle point 
in ip is located at 

TTl 

while the angles 791 and '§2 are not determined by the saddle point equations. The 
chiral condensate decreases linearly as a function of m until one reaches another second 
order phase transition for negative quark mass at the point where cos(p = —1. The 7 r° 
condensate can be calculated by switching on a small twisted mass as source term and 
setting it to zero after differentiation and taking the thermodynamic limit. Its value is 
non-zero in the Aoki phase. Thus parity and flavor symmetry are spontaneously broken. 
For Pi = (u = 0 the partition function is isospin symmetric. A source term corresponding 
to a charged pion condensate i± —)■ 0 yields a hysteresis and polarizes this condensate 
due to spontaneous symmetry breaking. The phase diagram at pi = ca = 0 is shown in 
Fig.gi). 

3.1.3. Non-zero Lattice Spacing and non-zero Chemical Potential. For imaginary 
isospin chemical potential and vanishing twist the phase diagram looks almost the same 
as the one at pi = 0, cf. Figs, [^i) and ii). The reason for this behavior is that 
sin'll = 0 minimizes the saddle point action for \rh\ < 16af. This hxes the charged 
pion condensates and their corresponding covariances at zero while the 7 r° condensate 
remains non-zero even in the absence of a source term. The flavor symmetry is broken 
explicitly by the chemical potential so that the would-be Goldstone bosons are massive 
in this phase which is labeled by The subscript “±” refers to the sign of the 

quark mass since the non-vanishing 7r°-condensate switches sign when m crosses the 
origin. Because this behavior is reminiscent to the Sharpe Singleton scenario (denoted 
by S.S.^) where the mass dependent chiral condensate jumps we refer to it as a Sharpe- 
Singleton-like scenario and abbreviate it by S.S.^. Note that both parts, and 

refer to the same phase although they are connected by a discontinuity along 
the m = 0 axis. This part of the phase diagram is shown in Figs. I^ii), iv) and v). The 
details are summarized in Tables [T] and [H 

At Pi = 0 we are in the Aoki phase with two exactly massless Goldstone bosons 
if |m| < 16a^. At a non-zero real isospin chemical potential in the related regime 
\fh\ < 16a^-|-pj/2 we immediately enter a Bose condensed phase with an isospin density 
that initially increases linearly with pi. The reason is that an exactly massless Goldstone 
boson with non-vanishing isospin charge exists in this phase. This phase is denoted by 
I. Note that for non-zero imaginary chemical potential, no Bose condensation takes 
place (phase II±), but the pions become massive because of the symmetry breaking 
by the chemical potential. The free energy in the phase I is quadratic in pi for small 
chemical potential so that we have a second order phase transition for \rh\ < 16af at 
pf = 2|m| — 32a^. 

In the regime \m\ > 16a^ -|-pi/2 we hnd the normal phase where the mass 

dependent chiral condensate S(m) is a constant (cos(p = 1 ) and all the other order 
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(i) 


gj<0 



sis." 
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oj 

ip 
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Figure 3. Phase diagram in the rh — a^-plane (Figs, (i) and (ii)) and in the gf-fh- 
plane (Figs, (iii) and (iv)) for 0 < w < 7r/2. In particular we choose (i) gf < 0, (ii) 
hi > 0, (iii) < 0 and (iv) > 0. The phase / only exists for real isospin chemical 
potential. For gf < 0 and < 0, the chiral condensate jumps when the mass crosses 
the origin and while the 7r°-condensate jumps for > 0. This case is related to the 
Sharpe-Singleton scenario. The curves of the second order phase transitions at gf > 0 
in hgures (ii) extend to infinity at m = ±/Ij/2sina;. There is no phase transition 
beyond these masses in the leading order expansion of the chiral Lagrangian. A finite, 
imaginary effective lattice spacing < 0 shifts the bifurcation of the phase transition 
curves away from the origin (see Fig. (iii)). 


parameters vanish. In contrast, for \m\ < 16a^ + fif/2 (phase I) we have 

I cos 9 ?I < 1 and cos-di = 0, (3-12) 

so that the 7r° condensate always vanishes and parity is not spontaneonsly broken. The 
flavor symmetry is broken explicitly by the chemical potential thongh. In the phase 
I the chiral condensate drops linearly in the qnark mass when approaching m —)■ 0 as 
in the Aoki-phase. The phase bonndaries are shifted down by the valne —/if/32, in 
contrast to the case /if < 0 cf. Figs, [^ii) and iii). In the presence of a charged pion 
sonrce the corresponding condensate is nonvanishing in the phase I. This part of the 
phase diagram is shown in Figs, [^iii), iv) and v). 
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At imaginary isospin chemical potential the saddle point equations dictate cos'll = 
—signm and cos(p = signfh cos ^ {fp G [0,7r/2]) so that the phases are determined by 
the minimum of 


VC 


01 cos 1 ^ 1 =—sign m 


—2|m| cos{ip — ca) + 16a^ cos^ p. 


(3.13) 


In Fig. we plot the solution of this equation for various values of a; G [0,7r/2] as a 
function of |m|/16a^. 

The Sharpe-Singleton(-hke) scenario extends to all of as well as masses since the 
twisted mass wipes out the second order phase transition to the normal phase. In fact 
the Aoki phase as well as is absent for 0 <a;<7 r/2as can be seen from the 

a; —)■ 0 limit of the curves in Fig. [Tj For > 0 the 7r° condensate changes from S to — S 
when the quark mass traverses zero while the mass dependent chiral condensate S(m) 
jumps for cP < 0. The two minima are separated by a potential barrier so that we have a 
hrst order phase transition when the mass crosses zero. The corresponding two parts of 
the phase are denoted by II± (again ± refers to signm = ±1). In this phase the chiral 
condensate is given by cosp and the condensate by sin(y9 so that the combination 
YP{m) + = TP is constant, cf. Fig. |^i) and Table The chiral condensate and the 

7 r° condensate only depend on the angle ca and on the ratio m/cP such that we have 
some kind of a “Silver-blaze property” in the isospin chemical potential. 

The situation changes when the isospin chemical potential is real. For 2|m| sinca < 
pf the effective potential also has a minimum at 


cos'll = 


msmcj 


(3.14) 


4/ij sin p ’ 

which yields the phase I. It shows up above the value of > —'fli/32 and between the 
masses — fif/2 sin u < fh < /ij/2sina;. The region where cP < —/if/32 is reminiscent of 
the Sharpe-Singleton scenario since the chiral condensate is jumping at m = 0 while the 
7 r° condensate does not because sin p vanishes for aP < 0. The two second order phase 
transition lines asymptote to m = Pfif/2 sin uj. Hence we do not have the situation as 
for /ii = a; = 0 where one can increase the effective lattice spacing, and independent of 
the value of the quark mass, one always enters the Aoki phase. For sufficiently large 
quark masses the system will always stay in the phase II± and will never enter the 
phase / if w > 0, cf. Fig. [^ii). However when increasing the isospin chemical potential 
fii the system will eventually enter the phase /, see Fig. [^iv). We underline that the 
described behavior can be caused by the fact that we only took the leading order of the 
chiral Lagrangian into account. 


3.3. Maximal Twist 

At real isospin chemical potential the cup-like structure of the phase transition lines 
becomes more and more rectangular shaped when increasing the twist angle. Eventually 
the smooth curves describing the phase transition lines develop a kink in the limit of 
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Figure 4. Phase diagram at maximal twist (w = 7r/2) in the m — a^-plane (Figs, (i) 
and (ii)) and in the gf — m-plane (Figs, (iii) and (iv)). We draw the phase diagram 
at imaginary isospin chemical potential ((i) /ij < 0), at real isospin chemical potential 
((ii) /ij > 0), at imaginary effective lattice spacing ((iii) < 0), and at real effective 
lattice spacing ((iv) > 0). Phase / only appears if hi > 0. The parts and 

jj^-tt/ 2 ) p];iase are separated by a Sharpe-Singleton transition (chiral 

condensate jumps when m crosses the origin) and extends to the phase 
where the 7r° condensate jumps. The case of real effective lattice spacing ((iii) P > 0) 
and maximal twist is the analogue of the case at imaginary effective lattice spacing 
and no twist, cf. Fig. [^(iv) and vice versa. Comparing the phase diagrams at no and 
maximal twist we observe the symmetry (w, to, , a^) —>■ (7r/2 — uj, fh, hi + 32a^, —ah) 
of the static Lagrangian, cf. the symmetry transformation 1^ . 


maximal twist {u = 7r/2), see Fig. |^ii) and iii). At the corners a new second order 
phase transition to the “normal” phase appears, which corresponds to the 

solntion (p = 7r/2 that develops in the limit a; —?• 7r/2 for < 0. 

In the regime |m| < 16a^ we find another phase Figs. 4 ii) and iii). In 

both kinds of phases we have a non-vanishing 7r° condensate. However only in 
this condensate becomes independent of the qnark mass and is some kind of a “Silver- 
Blaze property” in all parameters |/iip, \m\, and |ap while in have a linear 

dependence on m/16a^. Fnrthermore, the mass dependent chiral condensate S(m) is also 
non-zero in the phase becomes ±S at very small qnark masses exhibiting 

the original Sharpe-Singleton scenario [^. This is not the case in the phase 
If the effective lattice spacing of /V satishes of > —/ij/4, one does not cross the phase 
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Figure 5. The three pion masses at vanishing effective lattice spacing (a = 0) and at 
real isospin chemical potential {gf > 0) as a function of |to| (left plot) and of |)2i| (right 
plot). Notice that the squared pion mass is plotted in the right figure to underline the 
mass dependence. In this figure we show the functional behavior of the pion masses 
at zero chemical potential (blue curve in right figure) and at zero mass as a function 
of the chemical potential (blue curves in left figure). The neutral pion mass remains 
massless as a function of the chemical potential for zero quark mass. This diagram 


was already derived in 28,29 


jj^-tt/ 2 ) varying the quark mass but the phase / instead. 

Both phases and carry over to imaginary isospin chemical 

potential while the phase I does not exist in this region. The latter behavior was 
already the case for vanishing and hnite twist. Note that the whole phase diagram at 
maximal twist is a flipped version of the phase diagram at vanishing twist reflecting the 
symmetry (2.20), cf. Figs. and The Sharpe-Singleton transition lines inside the 
phases and along the m = 0 axis can be also found at vanishing 

twist. In the region > 0 and pf < 0 the condensate flips its sign when crossing 
with the quark mass the origin. At imaginary effective lattice spacing < 0 the 7r° 
condensate gets a kink at a second order phase transition and then drops linearly off 
while the mass dependent chiral condensate S(m) discontinuously switches the sign with 
the quark mass. The phase diagram is shown in Fig. and the details are summarized 
in Tables [T] and |2l 

The special case fif = —32af > 0 is related to the Aoki phase by the symmetry 
transformation (2.20). Then the parameters of Uq are only constrained by 

‘2iTTi 

sin cos'll = — (3.15) 

hi 


Hence we 

■l2 


have two massless Goldstone bosons when taking the limit flf = 
+ e). Note that this limit is sometimes different from the limit pf = 

This is the crucial difference between the original Aoki 

phase at a; = 0 and the one identihed at u = ti/2. 


hm(32|a 

e\0 

hm(32|a|^ — e), see section 

e\0 
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The masses of the pions are obtained by expanding the chiral Lagrangian to second 


order in the pion helds, as it is done in £ 2 , cf. (2.19). Thereby we only consider a real 
isospin chemical potential. 

The pole masses are defined by the energies B when the qnadratic form containing 
the pion fields, 


^k-ldki(B,p)ni 


(4.1) 


see Eq. (B.ll), where we have nsed Minkowski metric, becomes singnlar. In particnlar 


we have to determine B for which 
det II(B,p = 0) = 0. 


28 


(4.2) 


This calculation is worked ont in Appendix B The resnlts at vanishing and maximal 
twist are snmmarized in Table [H 

Physically it is clear that for a = 0 and snfficiently small pi (snch that the system 
is in one of the phases II± and III±) the chemical potential dependence of the pion 
masses is given by 

= 0) + g/ii, (4.3) 

where q is the isospin charge of the pions, compare with the linear dependence of the 
cnrves in Fig. The pions condense for pi > m^(/ii = 0) resnlting in one exactly 
massless Goldstone boson (phase I), since one of the angles in Uq is not determined by 
the saddle point eqnations. 

We 


At vanishing twist and lattice spacing the pion masses are well known 28,29 


recall the detailed dependence of the pion masses on the qnark mass and the iso spin 
chemical potential in Fig. [h In the “normal 


phase the mass relation (4.3) 


is exact with VB‘^m^q{iii = 0 ) = 2 |m| while in the pion condensed phase we have one 
massless charged pion and two massive pions which are the other charged pion and the 
nentral one. 

The next simplest case is = a; = 0 and > 0. This case was analyzed in 10,11 


In the Aoki phase, m < 16a^, the saddle point eqnations only fix the angle (p, with a 
nontrivial dependence of Uq on "di and '^2 remaining. This is consistent with having two 
massless pions in this phase at zero twist. At the phase transition point also the third 
pion becomes massless (see Fig. [^for pf = 0). At a finite twisting angle a; 7 ^ 0 all three 
pions become massive when the qnarks are massive (m 7 ^ 0) . The reason for this is 
that we never enter the phase / for a; 7 ^ 0 , /Ij = 0 and > 0 . 

When the Aoki phase is pertnrbed by a real isospin chemical potential, one of the 
charged pions condenses. Despite the symmetry breaking by the chemical potential, one 
of the two massless pions therefore remains massless in the phase I for all valnes of the 
twist angle, cf. Fig. [^c). When we are ontside the Aoki phase for pi = 0, there are no 
massless pions at pi = 0. However one massless pion always appears for snfficiently large 
isospin chemical potential when entering the phase I (see Fig. id). Also for < 0 we 
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\M 

^[32\W 

(d) 


IM 


Figure 6. The behavior of all three pion masses drawn in the m 7 r-|/i|-plane for various 
twisting angles uj and real isospin chemical potential {pf >0). In this figure we show 
the situations: (a) 0 < |m| < —16a^; (b) |m| > —16a^ > 0; (c) 0 < |m| < 16a^; and 
(d) |m| > 16a^ > 0. The vertical light green and dark green dashed lines indicate the 
positions of the phase transitions between the phases / and III± at no and maximal 
twist, respectively. The jump of the masses in figure (a) for uj = 7r/2 reflects the first 
order phase transition between the phases I and Moreover there is no phase 

transition in figure (c) for a; = 0. The horizontal marks on the y-axis for uj = 0,7^12 
are for orientation in the diagrams. In particular it should underline the generic form 
of these pion mass curves. 


find a transition to a pion condensed phase for pi > = 0), see Figs. |^a) and b), 

and[^ 

At maximum twist we enter an Aoki like phase at /ij = —32a^. Depending on 
how we approach this surface we have one or two massless pions, see Fig. |^a). This 
discrepancy reflects the fact that the transition from to the pion condensed phase 

/ is first order. This is shown in Figs, [^d) and f) where the blue curves correspond to 
Ut = lim(—32a^ — e) while the green ones are the one for uf = hm(—32a^ + e). The 

^ £->0 / o ri 

latter limit yields two massless Goldstone bosons which again reflect the fact that two 
angles in the unitary matrix Uq are not fixed, see also the Table 


4.1. Phase I (AfP > |m|sina; > Oj 

In this phase both p and "di are neither zero nor 7r/2. Then the U (1) symmetry of 
the microscopic theory is spontaneously broken resulting in a massless Goldstone boson 
corresponding to the angle ’§ 2 - The saddle point solution Uq is a combination of all three 
classical pion fields (vr*^ and tt^) because in the second order Lagrangian (2.19) all three 




























































Phase Diagram of Twisted Mass Wilson Fermions 


25 



VF>^.3 

32 W 


(C) 



Figure 7. Pion masses at real effective lattice spacing (fp > 0) in the TO^-|m|-plane 
for various twisting angles u. Shown is the generic behavior at vanishing isospin 
chemical potential (blue, solid curves) and at finite real isospin chemical potential 
(red, dashed curves). The vertical dashed lines indicate the positions of the phase 
transitions between the phases I and III± at zero (dark green) and maximal twist 
(light green). The transition between the phases I and II± at finite twisting angle is 
always between those two points. The kinks of the masses below these positions result 
from exact crossings of the masses and the ordering < m 7 r ,2 < w,r, 3 - The phase 
diagram at = 0 was discussed in [T^. 


modes are coupled. This results in a somewhat complicated expression for the masses 
mn,2 = \JCl - y/c2, fh^,3 = \JCl + \fc2 (4.4) 


with 

Cl= Pi + 

C2 = 

— 48a^/ij 


6m^ sin^ uj Qw? cos^ oj 


+ 


p\ (320^ + hi 

2 


+ 16a^ 1 - 


4m^ cos^ UJ 
(32a2 + /22)2 




1 - 


4m2 cos^ UJ 
(32a2 + filf 


9 ^ / 4m^ cos^ UJ AnP sin^ u 


+ 4 hT 


(32a2 + 


232 


Pi 


Idfh^cos'^o; sirP uj Arn'^ cos^ uj(jlf — iirP sirP uj) ', 


(32a2 + /22 


■2\A 


Pi 


yUj (32a2 + hj )2 


The third mass fh-n-^i = 0 vanishes. In Appendix B.l we show the details of the 
computation. 
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Table 3. Functional dependence of the three pion masses at zero (w = 0) and 
maximal (w = tt/2) twist. Note that we consider real isospin chemical potential 
(/tj > 0) such that the phase does not appear. At finite twist the expressions 

are quite complicated, cf. Eqs. (4.4) and (B.20). Hence we omit them in this list. We 
underline that we can approach the Aoki-like phase at maximal twist from two sides 
(gi = lim(—32a^ ± e)) yielding a jump in the pion masses due to the first order phase 

e \,0 


transition which happens at this phase. This cannot be done for the original Aoki 
phase at zero twist (uj = gi = 0) since gi has to be real and thus gf > 0. 
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0<-32a-</i; 
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32 |a|| 

0 < fij = -32 a- + 6 ' 

|m| 
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Figure 8. Pion masses at imaginary effective lattice spacing (a^ < 0) and at real 
isospin chemical potential (^j > 0) as a function of the quark mass |m| for various 
twisting angles w. Shown is the generic behavior for the following isospin chemical 
potentials: /tj = 0 (blue, solid curves in (a), (c), and (e)), /tj < —32a^ (red dashed 
curves in (a), (c), and (e)), fif = lim(—32a^ — e) (blue, solid curves in (b), (d), and (f)), 

/ij = lim(—32a^ + e) (green, solid curves in (b), (d), and (f)), and gf > —32a^ > 0 

e \,0 

(red dashed curves in (b), (d), and (f)). Note, the mass curve for is the same 
for gf = 32o^ ± e. The vertical grey, dashed lines indicate the positions of the phase 
transitions between the phases I and at no and maximal twist. Note that 

there is no phase transition for gf < 32a^ and uj -n 12. The kinks of the masses below 
the positions of the phase transitions for gf > —32a^ > 0 result from exact crossings 
of the pion masses and an ordering < mTr ,2 < w,r, 3 - The jump of the masses m,r ,2 
and m 7 r ,3 at g? = —32a^ reflects the first order phase transition between the phases / 
and . The curves for /ii = 0 were already derived in 

4 . 2 . Phase II± ( I cost^il = 1 and w 7^ 0, vr/2j 

In this phase the saddle point solution satisfies a reduced parameterization, 


19 


Uq = cos ip P i sin pr^. 


(4.6) 
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From the second order chiral Lagrangian (2.19) one can immediately read off that the 
7 r° mode decouples from the charged pions Moreover, the second order terms in vr® 
do not depend on the chemical potential fii. The corresponding masses are worked out 
explicitly in Appendix B.2[ From Eq. (B.20) we obtain (using that a = coslpo) 


2 2 |Sm| 32 a‘^C 2 ^ .2 2 ^ 

m^o = cos(v3o - tj) H-=^(sm cpo - cos (po), 


F3 




(4.7) 


where (po is determined by Eq. (A. 10). 


The chemical potential dependence of the charged pions is again given by 


relation (|4.3|). However the mass at vanishing isospin chemical potential may differ 

(see 


from the one of the neutral pion, i.e. (see Appendix B.2) 


2 I Xj?77/ I 

m^± (/q = 0) = ' 2 cos(^o - 


32a‘^C2 2 

cos (po- 


(4.8) 


■p2 ^ p2 

TT TT 

Additional details are worked out in [Appendix B.2 and can be read off Figs. and 

Table [H 


5. Summary 

Starting from the leading order chiral Lagrangian in the p-regime we studied the phase 
diagram of two degenerate twisted mass Wilson fermions as a function of the quark 
mass, the isospin chemical potential, the lattice spacing and the twist angle. This work 
extends previous studies of the phase diagram in the same parameter space but at 


zero lattice spacing and vanishing twist 28,29 or at zero isospin chemical potential and 


hnite twist and lattice spacing 19 . The phases are characterized by order parameters, in 


particular by the mass dependent chiral condensate and by the neutral and charged pion 
condensate, which enables us to distinguish the various phases. We consider both real 
and imaginary isospin chemical potential, and real and imaginary effective lattice spacing 
(because of the sign of the low-energy constants). Since the saddle point approximation 
can be analyzed for parameters in the e-domain we can always stay well away from the 


Roberge-Weiss transition 48 at imaginary isospin chemical potential. 


At zero chemical potential and at zero twist we can distinguish, depending on the 
values of the low-energy constants, two different possibilities in the approach to the 


chiral limit at non-zero lattice spacing. First, a transition to the Aoki phase 10 with 


a non-zero neutral pion condensate which spontaneously breaks parity, and second, the 


Sharpe-Singleton scenario 11 , where when the quark mass crosses zero, we jump to a 


different minimum separated by a potential barrier but with the same physics. 

A non-zero isospin chemical potential destroys the Aoki phase. Since the Aoki 
phase has two charged pions, a non-zero isospin chemical potential immediately leads 
to pion condensation so that the partition function becomes singular at /ii = 0. In the 
case of the Sharpe-Singleton scenario, the hrst order jump of the mass dependent chiral 
condensate persists at imaginary effective lattice spacing until the chemical potential 
reaches a critical value. Interestingly the Sharpe-Singleton scenario may extend to real 
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effective lattice spacing where not the chiral condensate bnt the vr® condensate may 
jnmp when the qnark mass crosses the origin. Moreover, we are always in the normal 
phase for large qnark mass. This phase is characterized by a non-zero valne of the mass 
dependent chiral condensate. At finite twist the same part of the chiral condensate is 
rotated into a 7r° condensate. 

At non-zero twisted mass, the phase diagram greatly simplifies and only two phases 
are possible. Phase / is connected to the pion condensed phase at zero lattice spacing 
when the chemical potential is less than the pion mass. The other phase II± is 
continnously connected to the “normal” phase. The two phases are separated by a 
second order phase transition. 

We have also shown that in the considered parameter space, the phase diagram 
has a symmetry which relates the mean field solntion at a given twist angle u to the 
complementary twist angle n/2—uj. We nnderline that this symmetry is not a symmetry 
of the mass spectrnm of the psendo-Goldstone bosons. In particular, this implies that the 
physical properties of the phase at maximum twist corresponding to the Aoki phase are 
different from those of the Aoki phase. Instead of two massless Goldstone bosons, we find 
either one or two massless Goldstone bosons at the critical value of the isospin chemical 
potential depending on how the phase boundary is approached, fif = lim(—32a^ — e) 

e\0 

and yhj = hm(—32a^-|-e), respectively. The pions in the phase / do not have good 

e\0 


isospin quantum numbers. Taking the isospin chemical potential away form its critical 
value does not result in pion condensation. 

The analytical results presented in this paper are useful for lattice studies of QGD at 
non-zero isospin chemical potential with Wilson fermions with zero or maximum twist. 
In particular they will be helpful to determine the region of the parameter domain that 
is connected to the physical point. Interesting extensions of this study would be to 
QGD with more than two flavors Is] , to lift the degeneracy in the quark mass as studied 


m 


50 52 for /ij = 0, and to QGD-like theories like two color QGD as it was recently 


analytically discussed by three of the authors 53,54 
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Appendix A. Derivation of the Phase Diagram 


In this Appendix we determine the solutions of the mean field equations. Different types 
of solutions are discussed separately in [Appendix A.l| and [Appendix A.2[ The case of 
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zero twist is worked out in Appendix A.3 and the case of maximum twist in Appendix 


In the thermodynamic limit we have to minimize Eq. (3.1) with respect to the 


variables ip G [0,7r] and iHi G [0,7r[. Before starting the calculation of the saddle point 
solutions let us discuss the Lagrangian in more detail. First, the Lagrangian does not 
depend on the third angle ^2 parameterizing the group SU(2) but this angle may be 
hxed by the source terms via spontaneous symmetry breaking. Second, the minima of 
the Lagrangian have to satisfy the following two conditions 

mcos<p>0 and mcos'di<0. (A.l) 

This hxes the signs of the trigonometric functions of the angles ip and "di which are 
sign cos ip = —sign cos'll = signm while sin(p = sin791 > 0 because G [0,7r]. For 


this choice the hrst two terms of Eq. (2.34) are negative which would be positive when 


choosing other signs. The other terms in Eq. (2.34) remain the same when switching 


the sign of the trigonometric functions. The third thing we want to point out is that 


the Lagrangian (3.1) is a quadratic function of the two parameters x = cos ip and 


y = sin ip cos 'di on the centered unit disk. Therefore we have either a minimum inside 
the disk < 1 equivalent to | cos'd!] < 1) or on the boundary = 1 

implying cos'd! = —signm). 


Appendix A.l. Inside the unit disk (\ cos'd!] < Ij 

The global extremum of the polynomial in the variables x = cos ip and y = sin ip cos'd! 

is 


2m cos oj 


X = 


and y = 


2m sin oo 




(A.2) 


32a2 + 

However it is only a minimum inside the unit disk x^ + y'^ < 1 (and thus a global 


minimum) if 


0 < 2]m] cosca < 32a^ + yf and 0 < 2]m] sinca < /ij. 


Additionally we have 

1 > \ll 


f 2 m cos 00 
\32a^ + fif 


> 


2]m] sinca 


hi 


2 — 


> 0 


(A.3) 


(A.4) 


which encodes the fact that we are inside the disc. We denote the phase corresponding 
to this situation by I. 

The angles take the values 

/n / 2 m cos u \ .(I) ( 2mcosoj 

ip^ ’ = arccos I ———^ I , 'd) = arcco 

y o2a ~\~ J 

such that the group element U G SU(2) is given by 


pf sin ipT) 


(A.5) 


= 


2m cos 00 
32a? + pj 




(A.6) 
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+ ^ 


2m sin oo 




4m2 cos^ oj Am? sin^ oj 


3 — 


Am? cos^ OJ Am? sin^ oj 


' 1 - 




(32a^ + 'il\y 

with ^2 G [0, 27r]. The free energy is then equal to 


(32a2 + ^2)2 
2m sin a; 

tI 


Ft 




/ij 2m? sin^ oj 2m? cos^ oj 


2 fil ?>2a? + fll 

Again we want to underline that the angle ^2 can be hxed by the source terms j±. 


(A.7) 


Appendix A.2. On the unit circle (\ cos'd!! — ^) ^ 7^ 0,7r/2 

The extrema on the boundary {x,y) = (cos</9, sin<y9) = (sign m cos sin are given by 
the equation 

16a^ cos ^ sin ^ — |m| sin ca cos ^ — |m| cos casing = 0 (A.8) 

which is the derivative of the chiral Lagrangian 

VCq = 16a^ cos^ ip — 2\m \ sin a; sin ^ — 2|m| cos a; cos (A.9) 

Note that this case is independent of the chemical potential such that the solution has 
a “Silver-blaze-property”. 

Let OJ 7^ 0,7r/2. Then the solution of the transcendental equation can be compactly 
written as an integral, 

|m| \ 


cosp = F 


Wa^J Jo ® 


mK^^l — j/2 cos a; — j/sinca) — 16a^y^/l — dy, (A. 10) 


where 0 is the Heaviside step function. One can easily show that Eq. (A.10) is the 


solution of Eq. (A.8) via integration by parts and taking into account that there is only 


one solution in the interval cos93 g] 0, 1[. Thus the angles freeze out at 

= arccos signm J 0 ^|m|(\/l — y‘^ cosca — ysino;) — 16a‘^y\/l — y^^ dy 
= arccos (—signm), 

and the corresponding unitary matrix is 


(A.ll) 


= signm J 0 ^|f«|(\/l — cos a; — j/sinca) — 16 a^y\/l — dy'l2 (A-12) 


— zsignmW 1 — 


Uo 


0 ^|m|(\/l — 2/2 coscj — ysino;) — lQa?y^/l--^'^ dy 


Ts- 


This solution corresponds to the phase labelled by “JJ±” where the subscript denotes 
the sign of the quark mass. The free energy becomes 

V = 8®^ cos2(p^^^^ — 2\m\ cos{(p^^^^ — ca) -|- 8a^. 


(A.13) 
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The complete region to find the phase II± is given by 


hi > 2 \m\ sincj > 0 and 


2 |m| sin a; 


> Wi- 


2 m cos u 


>0 or ( 2 |m| sin a; > hj) 


hi V V + hi ^ 

(A.14) 

Collecting everything we snmmarize that there are two phases, I and II±, at finite 
twist angle w 7 ^ 0 , 7 i/2. 

Appendix A.3. The limit of zero twist (co ^ 0) 

Let us now consider the limit of the two phases when setting the twist angle uj to zero. 
The first phase I will be in the region 

hi > 0 and 0 < 2 |m| < 32a^ + pj 

with the unitary matrix 

2 m 


(A.15) 


U^%=o = 


2 >2a? + h? 


1 

0 


(32S2 + /lf)2 


0 


and the free energy 


V h I — 

V n-o |a;=0 — 


Ml 


2 m? 


(A.16) 


(A.17) 


2 32 a ?+ -^1 

This phase meets the hi = 0 plane where the Aoki-phase is. 

What happens with the phase //±? In this case the solution of the transcendental 


equation (A. 8 ) breaks up into two branches, 

m 

smo 9 = 0 or cos<h= — 

loa"* 


(A. 18) 


The first branch corresponds to the region 

(hi > 0 and 2 |m| > 32a^ + fij) or (O > hj and \m\ > 16a^) , (A.19) 

and the second branch is found in the region 

hi < 0 and 16a^ > \m\ > 0 . (A. 20 ) 

We always have the solution tp = 0 which only dominates if we do not reach the solution 
(p = arccos(|m|/16a^) only occurring in the region (A.20). In both regions the correct 


solution is given by Eq. (A.IO). 


The first region ( |A.19 ) yields the unitary matrix 

= signmll2 (A. 21 ) 

and the free energy is given by 

= Wc? - 2\m\. (A. 22 ) 

Note that the phase denoted by (corresponding to signm = ±1) is continuously 

connected to the phase I implying a second order phase transition. 
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In the second region (A.20|) the unitary matrix becomes 

- *signm^ 


IGa^ 

and the free energy is given by 


/ rh 

I IGa^ 


T3 


VC 


iii)\ 

0 11^=0 


m 


16a^ 


(A.23) 


(A.24) 


At vanishing isospin chemical potential the phase meets the phase / at the Aoki 

phase which only exists in the (/ij = 0)-plane. Starting from a hnite twist angle the 
phase II± is continuously approached from the phase I. This is not true anymore for 
O' = 0 where it is a first order phase transition. To see this, one has to take the limit 


/ij —)■ 0 in Eqs. (A.16) and (A.23). 


We want to underline that signm has to be replaced by sign jo in Eq- (A.23) when 


we do not approach this phase from a hnite twist but at zero twist with the source term 
jo. Thus the order of the limits a; —?• 0 and jo —?■ 0 is crucial in this phase to determine 
the correct saddle point equation. 

The Aoki phase only exists at = 0 and is quite particular since the angle 'di also 
drops out while sin (f 0 such that the unitary matrix still depends on this angle 


with the free energy 


/l / m ^ 

2 

cosdi 

sin hi 


1 

sin hi 

— cos hi 


VC 


(Aoki) I 


m 


| t »;=0 


16 a 2 ' 


(A.25) 


(A.26) 


As a function of it connects the phases / and by a hrst order phase transition. 

Finite pion condensates can be created via the three source terms jk, k = 1, 2, 3. 


Appendix A.f. The limit of maximal twist (uj —)■ 7r/2j 
At maximal twisted mass the phase I yields a unitary matrix 

U ^%=./2 = 1 


and the free energy 

/ /ij 

in the region 

hi > 2|m| > 0 and 32a^ + pj > 0. 


2 m 




T- 


4^2 

hi 




T- 


4fh2 

hf 


2 m 


(A.27) 


(A.29) 
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In the phase II± the equation for ip again splits into two branches for u ^ 7 t/2 
(see Fig. [^. One branch with cosp = signm exists for \m\ > —16a^ and the other 
branch with sin(p = —|m|/16a^ is the minimum for 0 < |m| < —16a^. Taking the limit 
a; —)■ 7r/2 of the regions (A. 14) we thus obtain the following four regions 

{Pi < 0 and — 16a^ > \rh\ > O) or (—16a^ > \rh\ > 0 and — 32a^ > pi > O) (A.30) 
and 

(pi < 0 and \rh\ > —16a^) or (2|m| > pi > 0 and \rh\ > —16a^) . (A.31) 

Both possibilities are covered by the solution (A.10). 

When approaching the limit of maximal twist from a twist angle u < ti/2 in the 
phase II± we only smoothly reach the unitary matrix 

2 


= signm 

with the free energy 



m 




VC 


FT I 


\(jJ=7t/2 


m 


+ 16 a" 


(A.32) 


(A.33) 


in the region (A.30) if the quark mas s m r emains finite. However the other phase 
(denoted by JJ/^^A)) region (A.31) is approached by a second order phase 

transition. In this phase the unitary matrix is given by 

= -isignmrs, 

and the free energy is equal to 


VC 




|a;=7r/2 = -2|m|. 


(A.34) 
(A.35) 

Again the labelling reflects the sign of the mass signm = ±1. 

Similar to the case of vanishing twist we can find a different vacuum polarization 
in the phase interchanging the order of the limit u —?• 7r/2 and —?• 0. 

In the case that we do not regularize with a finite twist but with the source we have 
to replace signm by signmv in Eq. (A.32). 

Due to the symmetry (2.20) the Aoki phase at ca = 0 has a counterpart at ca = 7i/2 
which yields the saddle point 

= ^Z 


m 


16a2 

and the free energy 


D + 


/ m ^ 

2 

cos'll 

sin'di 

\16a?y 

} 

sin -di 

cos'll 


VC 


(Aoki) I 


|t»;=7r/2 


m 

1^ 


+ 16a" 


(A.36) 


(A.37) 


We already absorbed the phase shift resulting from the symmetry transformation (2.20) 
in the angle '^2 € [—tt, tt]. The phases / and j^gg^^ each other at this Aoki- 

like phase at pi = —32a^ > 0. The transition from / to jg eligcontinuous, cf. 

Eqs. (A.27) and (A.32), and is thus a first order phase transition. 

We find a finite charged pion condensate and chiral condensate in the Aoki-like 
phase at maximal twist when switching on the source terms j± and m^, respectively. 
Also a finite twist a; —)■ 7r/2 can align the chiral condensate. 
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To compute the pion masses we employ the vector notation of the Pauli matrices 


T = 


n 

T 2 

T3 


(B.l) 


Then any element in the Lie algebra su(2) can be written as 
A = AWj = (A, t) G su(2), with A= 

A3 


(B,2) 


where is the Euclidean scalar product in ^nd |.| is the modulus of the 

corresponding norm. Note that the norm itself will not be a map to non-negative real 
numbers if the vector is complex. Therefore we need the modulus. The Pauli matrices 
themselves are given by 

Tj = {ej,T) esu{2). (B.3) 

The chosen notation fulhlls the following rules 

AB = {A,t){B,t) = {A,B)l2 + z{Ax B,t), (B.4) 

trA = tr(A, r) = 0, (B.5) 

tiTjA = tr(ej, 7:)(A, r) = 2(6^, A) = 2 Aj, (B.6) 

for all A, B G su(2). 

With the help of this notation we decompose the matrix 

Uo = at2 + i{l3,T), = al2-t{l3,T), (B.7) 


where the components a, (5^ G [—1,1] with + \fd\^ = 1 can be read off from Eq. (2.31). 
The pion helds are 

n = WTj = (n,7:), = {W)*Tj = (n*,?:). (b.s) 


In this notation the six terms in the chiral Lagrangian C 2 , see Eq. (2.19) read 
tr nn^ 


trn[ntf/t,r 3 ]_f/o 


= 2|n|2, 

= — 4«(n* X n, + aft X §2 + {{3, ef)ft)), 


(B.9) 


tr[n,r3]_[nt, f/oVgt/o]- = 8(|^p - - 2 {P,e^y)\Il\^ + 8 {a^ - |^P)|(n, 63)^ 

+ 16a(n*,e3)(n,^ X 63 ) + 16(^,e3)(n*,e3)(n,^), 

tr(e*‘^'^^17o + = 4(acosa; — (/d, 63 ) sinci;)|n|^, 

|tr(t/o-f/d)n|2 =i6|(n,^)|2, 

tT{Uo + ?7d)nn'f tr(t/o + U^o) = lea^inp. 

The third term is generally complex and only becomes real after summing over all 
energies E and momenta p. Note that all terms are rotation invariant in the 6^-63 
plane so that we can rotate the vector (3 to ff = + /dgCg with = /df -f /d|. 
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The aim is to find the dispersion relations from the Lagrangian £ 2 , in particnlar 
we have to find the energies E = Ey/VF^^ > 0 for which the pion propagator diverges. 
Given the bilinear form in the pion fields 

1 


£2 = 


21 /F 2 


5 ] ntff(B,p)n, 


(B.IO) 


p,E >0 


we have to find the energies for which the determinant of 


HiE,p) = 


Pk^ 


E^ + pj (2/3_|_ — 1 ) + 2 m(a coscn — sincn) — 32a^a^ ) 11 


+ 


32a‘^^l 2ifiiE{l - f3l) 

—2i'/2iE{l — (df) 0 

((/ij + 32a^)/33 — 2ta'fLiE)f3± {—PiCt — 2i'jliE(d3)(dj 
= Ho + 2 ifiiEHi + 


E^)to 


((pi + 32a^)/?3 + 2tafiiE)/3± 
{-pfa + 2 ifiiEldo)ld± 
/f?(l-2/?i) + 32a2/?| 

(B.ll) 


vanishes. Using that the real vector —otfd±^ (1 — jdf)) is in the kernel of the 

antisymmetric matrix Hi the determinant of H can be rewritten as 

det H{E,p) = {pkf - E^f + tiHo{pkP^ - E^f - A\u\^plE\pkf - E^) 






{Pkf^ — E'^) — 4:{u, Hou)fLf E^ + det Hq. 


(B.12) 


The matrix Hq is the cnrvatnre of the potential and is therefore positive semi-definite 
at the saddle points. The norm of u is simply y/d — ftf while the other terms are qnite 
lengthy. 

Note in the phase II±, and the vector (dWe^, so that ld± = 0, 

Appendix B.2[ This results in a simplification of the Lagrangian £ 2 . For the phase 


see 


/, a simplification is that we always have an exact zero mode, see Appendix B.l 


Appendix B.l. Phased 

In the phase / we have (d± 7 ^ 0, and the parameters of Uq are given by 


2 m cos uj 


a = 


) 1^3 — 


2 m sin u 




1 — 




^l 


(B.13) 


32a2 + pj 

With these values we can simplify H to 
H{E,p)= (pk^-E^'^t, 

{fij + 32a‘^)/dl 2zpi£(l - 

- 2 zpi£(l - (df) ^ pfldl ^ 

{{pi + 32a^)(do — 2ta'jliE)(d± {—pia — 2ijliE(do)ld_ 

Then the determinant of Hq vanishes, i. e. det Hq = 0. Therefore one pion mode is 
massless = E^^^i^p = 0) = 0, see Eq. (B.ll The corresponding mode is given by 

(B.15) 


(B.14) 


+ 


((Pj 320^)jdo + 2iafiiE)/d_i 

{-fifa + 2ijliE(dii)jd±_ 

pi{l-Pl) + 32a^Pl 


n^,i cx -IdoTi + aT 2 + {l-a^ - IdiY^^To 


At the phase boundary the coefficient of vanishes and the mode 11^ 1 becomes a linear 
combination of the charged pion modes 11 + and n_ with a mixing angle determined 
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by the twist angle. At the phase boundary to phase II, this mode joins the mode that 
becomes massless starting from phase II. 

To find the other modes we have to solve the remaining quadratic equation which 
leads to the masses 


= 




(B.16) 


with 


Cl = ^(triJo + 4|M|^/ii), 

c, = - 4(«. 


(B.17) 


see Eq. (4.4). The corresponding modes involve quite complicated expressions which we 
omit. Nonetheless let us emphasize that they mix all three pion modes 7 r° and 
the phase boundary they become the vr^ mode and the massive charged pion. 


At 


Appendix B.2. Phases //±, and 


In the phase II± as well as have /3± = 0. Therefore the matrix (B.ll) 

simplihes to 

H= — hj + m(a cosw —/ds sino;) — 32a^a^') I 3 (B.18) 


0 21111E 0 

—2ifliE 0 0 

0 0 yUj + 32a^/3| 

From this expression we can simply read off the dispersion relations 


Et^o = \JpkjP + 2 m{a cos oj — sin oj) + 320 ^ (/d| — cn^), (B. 19 ) 

E^± = \/pk]P + 2 m{acosu; — /^ssino;) — 320^0,2 j- 

corresponding to the pions 7 r° = tsIIs and oc T=pn±, respectively. Recall the relations 
a = signm|a| and = —signfhVl — a^. Then the pion masses defined as the rest 
energy E{p = 0) = in.,, is 


m^o = 2|m|(|Q;| coscj + \/l — sinca) + 32ad{l — 2a‘^), (B.20) 

m^± = 2 |m|(|Q!| coscj + \/l — sinca) — 32a‘^a‘^ ± 'jli. 

At the phase boundary with the phase I one of the charged pions becomes massless and 
becomes degenerate with the massless pion from phase I. 

For the particular phases one has to plug the corresponding values for a into these 
expressions. Note again that only for pj > 0 the masses remain positive semi-definite. 
For yuf < 0 the masses become complex which is quite natural when considering the way 
the chemical potential affects the propagation of quarks. 
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